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Abstract 

The goal in this paper is to demonstrate a new method for 
constructing global-in-time approximate (asymptotic) solutions of 
(pseudodifferential) parabolic equations with a small parameter. 
We show that, in the leading term, such a solution can be con- 
structed by using characteristics, more precisely, by using solutions 
of the corresponding Hamiltonian system and without using any 
integral representation. For completeness, we also briefly describe 
the well-known scheme developed by V. P. Maslov for constructing 
global-in-time solutions. 



Introduction 

The goal of the present paper is to present a new approach to the con- 
struction of singular (i.e., containing the Dirac 5- function as a summand) 
solutions to the continuity equation and to show how these solutions can 
be used to construct the global in time solution of the Cauchy problem 
for Kolmogorov-Feller-type equations with diffusion, potential and jump 
terms. It is well known that the asymptotic solutions of the Cauchy 
problem for linear equations with a small parameter e > can be con- 
structed by the WKB method [27]. In the framework of this method, 
the initial partial differential equation is reduced to a system of equa- 
tions consisting of the Hamilton-Jacobi equation and several transport 
equations. All these equations can be solved under the assumption that 
the Hamilton system has smooth solutions corresponding to the above- 
mentioned Hamilton-Jacobi equation (the trajectories of the Hamilton 
system fiber the phase space). In general, this ensures only the existence 
of the classical solution in small with respect to time. In this case, if, for 
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example, the Hamilton function is time-independent, then, on the time 
intervals where the Hamilton-Jacobi equation has a smooth solution, 
the Cauchy problem for this equations (as well as for the corresponding 
transport equations) is invertible in time. 

It is well known that the above system of equations (the Hamilton- 
Jacobi equation and several transport equations) arises in the construc- 
tion of WKB solutions of the form Uas = exp{j^S{x,t)}Lp{x,t) for wave 
equations and Schrodinger-type equations [27] and in the construction 
of approximate solutions of the form = exjp{—^S{x,t)}ip{x,t) for 
parabolic equations [23]. 

We note that, in the first case (for stationary symbols), the property 
of being invertible in time is typical not only for solutions of this type but 
also for more general solutions because of the properties of the equation 
itself; but in the second case, the general solutions of the Cauchy prob- 
lem do not have this property, and the fact that smooth "WKB-type" 
solutions ex.-p{—^S{x,t)}(p{x,t) are invertible in time distinguishes this 
class of solutions from the other solutions. 

Now wc explane the time-invertibility condition in more detail and 
describe (below in this paper) a method for constructing such solutions. 

First, we note that the function S{x, t) can be defined as the point- 
wise limit S{x,t) = —lim^^Qln Us{x,t), where is a solution of the 
equation of parabolic type with a small parameter [24]. 

If the function S{x, t) thus defined exists and is smooth, then it is a 
solution (classical) of the Hamilton-Jacobi equation [24, 25] 



There is a well-known exact formula expressing the solutions of this 
equation in terms of the trajectories of the corresponding Hamilton sys- 
tem 



The fact that the solution of Eq. (0.1) is smooth for t = to means 
that the Lagrangian manifold obtained by a shift of the initial man- 
ifold Aq = {x, If |t=o} along the trajectories of (0.2) can be uniquely 
projected on i?" for all t £ [0,to]. 

Since the shift along the trajectories of (0.2) is invertible in time, 
this implies that the resolving operator of Eq. (0.1) is also invertible in 
time. Under the same conditions of unique projection, the solution of 
the transport equation for the function (p{x,t) (the amplitude) is given 
by the formula 



St + H{S^,x) = 0. 



(0.1) 



X = H{p, x) 



P = -H^{p,x). 



(0.2) 



ip{x,t) = 



C 



(0.3) 



VJ(M)' 
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where C is a constant along the projection of the trajectories of (0.2) on 
and J{x, t) is the Jacobian of the mapping of shift along these projec- 
tions (the uniqueness of the projection mapping implies the uniqueness 
and invertibility of the shift mapping). 

Under the above conditions, formulas (0.3) are also invertible. Ev- 
erything said above can be illustrated by the following simple example. 
Let us consider the heat conduction equation 



= e-^°(^)/^</?o, (0.4) 

t=o 



where the function So{x) > is assumed to be smooth and bounded 
together with its derivatives and (p G Cg", for example: 



So{x) = Hi 
Jo 



+ tanh z)dz = x + In coshx. 



Here the following addition condition is satisfied: 

The Hamilton function corresponding to (0.3) is 

H{p, x) = 

and, respectively, the trajectories of system (0.2) have the form 

a;(a;o, t) = xq + 2tp{xo, t), p{xo, t) = p{xo, 0) = = 1 + tanh xq. 

UXq 

(0.6) 

Obviously, condition (0.5) implies the unique globally-in-t solvability 
of the equation 

. dSo 
X = xo + 2t- — 
oxo 

in ixo and hence the global solvability of the corresponding Hamilton- 
Jacobi equation and the transport equation in the class of smooth func- 
tions. 

Further, for Eq. (0.3) it is possible to write the Green function, which 
is a solution of Eq. (0.4), 

such that G{x,^,0) = 6{x — It is well known that one can write 

u = J G{x,^,t)u{^,0)d^. 
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After the change t —t, we obtain the Green function for the inverse 
heat conduction equation, and it is easy to verify that 

u\t=o + O(e^) = J Gix, e, -t)e-^«'*)/^^(e, t) d^, (0.7) 

where S{x, t) and ip{x, t) are solutions of the Hamilton-Jacobi equations 
and transport equation 

^ + 2—^ + — -^ -0 
dt dx dc dx^^ 

at time t satisfying the initial conditions 

S\t=o = So{x), (p\t=o = ipo{x), 
where A/" > is an arbitrary positive number and 0{£^) is such that 

Equality (0.7) can be verified using the Laplace method. The last 
equality means that the estimation 0{£^) is not quit suitable here. 

There is an absolutely different situation if (0.5) is replaced by the 
inequality 

at least for a certain value of xq, for example, 



/•X 

So = (1 — tanh z)dz = a; — In coshx. 
Jo 

In this case, the jacobian of the mapping of shift along the trajecto- 
ries X = xo + is zero at t* = (min |2^^|)~''^ (in the example, we 



^2 



have = — cosh ^(xq) and t* = 1/2), where the maximum is taken 



dx, 







over all points at which the inequality (0.8) is satisfied. In our example, 
this is Xo = 0. For t > t*, the Lagrangian manifold shifted along the 
trajectories (0.6), 

AI = {x = Xq — 2i(l — tanh(xo)), p = 1 — tanh(xo)}, 

forms an ^-shaped curve, and there three points of this curve over a 
certain point x in the plane {x,p). These three points are associated 
with three (local) "WKB-type" solutions 

uj = e-^^^'''*y'^j{x,t). 
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It is clear that their hnear combination 

3 

is also a solution, i.e., satisfies the equation with the same accuracy as 
each of the functions Uj, j = 1,2,3. But the functions themselves are 
not equivalent. 

For example, it is clear that the inequality 

Si{x,t) > S2{x,t) 

holds at a certain point x, then the "WKB" solutions ui and U2 at the 
point X satisfy the relation 

= U20{e^), (0.10) 

where iV > is an arbitrary number. This follows from the fact that 
the difference (^i — S2)\x in parentheses in the exponent is positive. 

Thus, at each point in formula (0.9), it is necessary to choose the term 
where the function Sj is minimal. Such a choice leads to an expression 
of the form 

tx = e-<^(^'*)/^(^(a;,t), (0.11) 

where = (t>{x,t) = miax{Sj{x,t)}. It is clear that expression (0.11) is 
the leading term of the approximate solution. But its substitution into 
the equation in order to verify that it is a solution in a certain sense 
is not a trivial problem, because the function ^(x, t) thus determined is 
not smooth but is only continuous with a bounded first derivative. 

Of course, we can avoid this difficulty if we first calculate terms in 
(0.9) and then pass to (0.11). Such a construction, which takes into 
account the fact that the functions Sj and ipj can loose smoothness at 
the points where the Jacobian of projection mapping of the Lagrangian 
manifold on i?" is zero, was proposed by V. P. Maslov [26] and is ideolog- 
ically similar to the construction of the Maslov canonical operator (the 
Fourier integral operators) [8, 10]. A version of Maslov's construction 
was proposed in [28, 9]. 

The above procedure is quite well for constructing an asymptotic 
solution of the Cauchy problem; some estimates of the discrepancy be- 
tween the asymptotic and exact solutions were also obtained. 
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But attempts to apply these constructions to solve the inverse prob- 
lem meet insurmountable difficulties, and the details of the construction 
based on the use of integral representations do not play any role here. 

The problem consists in information that we need to obtain from the 
solution at t = T and the use to reconstruct the solution for < t < T. 

Indeed, as was shown above, if there are singularities in the pro- 
jection of the Lagrangian manifold on then the solution can have 
the form (0.9) at certain points, but as was shown above, some of the 
terms in this sum are "infinitely" small compared with the other. We 
can "measure" only expressions of the form (0.11) which do not contain 
information about the parts of the Lagrangian manifold corresponding 
to the "infinitely" small terms of the solution. But if we move backwards 
in time, then these (unknown) parts of the Lagrangian manifold can get 
into the domains that is uniquely projected on i?", and then they are 
responsible, in the projection, for the principal part of the solution. 

Thus, 

(i) If for t = T the function (/){x,t) = lim£^o(~^lii''^e) is smooth, then 
the principal part of the solution of the Cauchy problem can be 
reconstructed for < i < T in the entire space R^. 

(ii) If for t = T the derivatives of the function (p{x, t) have singulari- 
ties (discontinuities), then, for < t < T, the principal part of the 
function Ug (the solution of the Cauchy problem) cannot be recon- 
structed in the entire domain, because there are domains where it 
is impossible to reconstruct the solution. 

In our future considerations for non smooth case we will use the 
relation between transport and continuity equations. This relation be- 
tween the solutions of the continuity equation and the system consisting 
of the Hamilton-Jacobi equation plus the transport equation has been 
well studied before in the case of a smooth action functional. 

Let the velocity field u be determined as the family of velocities of 
points on the projections of the trajectories of the Hamiltonian system 
corresponding to the Hamilton-Jacobi equation. In this velocity field, 
as it was mentioned by E. Madclung [13], the squared solution of the 
transport equation satisfies the continuity equation 

Pt + iV,up) + ap = (0.12) 

with some additional term ap, which is defined below (a is equal to 
if the Hamiltonian is formally self-adjoint). The main obstacle to the 
extension of this correspondence globally in time is the fact that in 
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general the solution of the Hamilton-Jacobi equation are smooth only 
locally in time. The loss of smoothness is equivalent to the appearance 
of singularities of the velocity field mentioned above. Till the recent time 
there was no method for constructing formulas for solutions of continuity 
equation for a discontinuous velocity field. In Madelung's approach the 
divergent form of the continuity equation (in difference with transport 
equation) is very important property allows precisely to introduce a 
concept of global solution in spite of singularities in the velocity field. 

In the present work we generalize Madelung's approach to the case 
in which the singular support of the velocity field is a stratified manifold 
transversal to the velocity field trajectories. This holds, for example, in 
the one-dimensional case under the condition that, for any t G [0, T], 
the singular support is a discrete set without limit points. 

The class of solutions constructed in this way admits the motion 
forward and backward in time. Here we discuss only the construction 
as itself, the invertibility problem we plan to discuss in our next paper. 

1 Generalized solutions 

of the continuity equation 

Here we follow the approach developed in [1], where the solution is un- 
derstood in the sense of an integral identity, which, in turn, follows from 
the fact that relation (0.12) can be understood in the sense of distribu- 
tional space P(M"^^). The first step in this way has been done in [14], 
see also [15], [16], where the approach based on smooth approximations 
of the solutions was used. 

We specially note that the integral identities in [1] can be derived 
without using the construction of nonconservative products [2, 4] of the 
nonsmooth and generalized functions (or measure solutions [5] ) , and the 
value of the velocity on the discontinuity lines (surfaces) is not given 
a priori but is calculated. In the case considered in [1], the integral 
identities exactly coincide in form with the identities derived using the 
construction of a nonconservative product (measure solutions) in the 
situation described at the end of above introduction, which we shall 
now make more precise. 

First, we consider an n — 1-dimensional surface 74 moving in R", 
which is determined by the equation 

It = {x;t = t/j{x)}, 

where tp G C^(M"), and VV' 7^ in the domain in where we work. 
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This is equivalent to determining a surface by an equation of the 
form 

S{x,t) = 

(5 G in both variables, S{x,t) = 0, 'Vx,tS\s=o 7^ 0) under the condi- 
tion that 

dS , 

We remain that the situation with ^ = 0, can also be covered 
by making the change of variables x'^ = Xi — Cit with appropriately 
chosen i = l...n, solving the problem with the moving surface and 
then returning to the original variables. Possible generalizations are 
considered later in this section. 

Next, we assume that C belongs to C^{M."- x R\.). Then, by defini- 
tion, 

{d{t - V'(x)), C(x, t)) = [ C{x, ^{x)) dx, 

where 6 is the Dirac delta function and (,) is the distributional pairing 
(with respect to the variable t G and x G M"). 

Let S{t — tp{x)) be applied to the test function 77 G Co°(M"), then 



{6{t-il;{x)),r]{x)) = r)du^, 
J It 



where du) is the Leray form [6] on the surface {t = ip{x)} such that 
—dipdio^ = dxi . . . dxn- 

One can show that (see [1], [6]) 



First, we assume that the solution p to Eq. (0.12) has the form 

p{x, t) = R{x, t) + e{x)6{t - ip{x)), (1.1) 

where R{x,t) is a piecewise smooth function with possible discontinuity 

at {t = tp{x)}: 

R = Ro{x,t) + H{t - 'tp{x))Ri{x,t), 

e G C(IR") and has a compact support, i/j G and VV' / for x G 
suppe, and H{z) is the Heaviside function. 
It is clear that the term 

e{x)S'(t-4^{x)) 
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appears in (0.12) if we differentiate the distribution 6{t — '^{x)) with 
respect to t. Hence it is necessary to have in (0.12) 

(V, pu) = —e{x)5'{t — ip) + smoother summands, 

since V6{t — ^) = —VipS^t — tp). Then we must have 

eVip 

pu = -—-—7jS(t — ijj) + smoother summands. 

I V'f/'l^ 

Now we formulate an integral identity, defining a generalized solution 
to the continuity equation. 

We set r = {{x,t);t = this is an n-dimensional surface in 

M" X M]j.. Let 

u{x, t) = uo{x, t) + H{t — ijj)ui{x, t), 

where ijj is the same function as before, and U(),ui £ C(M" x M^). 

Let us consider Eq. (0.12) in the sense of distributions. For all 
C{x,t) e C^{W X R\), C(x,0) = 0, we have 

+ (V,H,C) = -{p,Ct) - {pu,VO- 

Substituting the singular terms for p and pu calculated above, we come 
to the following definition. 

Definition 1.1 A function p{x,t) determined by relation (1.1) is 
called a generalized 5-shock wave type solution to (0.12) on the surface 
{t = ip{x)} if the integral identity holds 



{RCt + {uR, VC) + aRQ dx dt 



for all test functions C(x,i) G P(M"xR^), C(x,0) = 0, = (||||,V) + 
We have also the relation 



/ 



C{x,^)dx = l^^J^C{x,t)dx. 



iRn IVV'I dn 



We note that the vector n± is orthogonal to the vector {Vip,—1), 
which is the normal to the surface F, i.e., lies in the plane tangent 
to F. 
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We can give a geometric definition of the field The trajectories 
of this vector field are curves lying on the surface F, and they are or- 
thogonal to all sections of this surface produced by the planes t = const. 
Furthermore, it is clear that the expression is an absolute value of 
the normal velocity of a point on 7^, i.e., on the cross-section of F by 

the plane t = const, and the expression • =^ Vn is the vector of 
normal velocity of a point on 7^. Thus, we have another representation: 



|V^| dn± 



Cix,t)dx = J^e(^iVn,V) + —y{x,t)dx, 



where Vn = 7r*(f„), v„ is the normal velocity of a point on 7^, and tt* is 
induced by the projection mapping tt : F — )• i?". 

It follows from the latter definition that the following relations must 
hold: 

Rt + (V, Ru) + aRC = 0, for ah points {x, t) ^ F, 

([i?]-|VV|[i?Un])+ (^^) r^"°' for aU points {x,t)eT, 

(1.3) 

The last relation can be rewritten in the form 

)CE + ^E= [Run] \Vi;\-[R], (1.4) 

where E = e/|VV'|, the factor /C = (V, |^^) = divi^ {v is the normal on 
the surface {t = ip{x)}) and, as is known, is the mean curvature of the 
cross-section of the surface F by the plane t = const, ^ = (j^^, V). 
Now we assume that there are two surfaces 

ri = {{x,ty,t = iiJi{x)} 

in R" X M^_, i = 1, 2, whose intersection is a smooth surface 

7 = {(ar,t);(i = Vi)n(t = V2)} 

belonging to the third surface F(-3) = {{x,t);t = ip3{x)}. Further, we 
assume that the surface F(3) is a continuation of the surfaces F^*) in 

(i) 

the following sense. We let n\ denote the curves on the surfaces Fj 



and we assign that each point (x,f) on the surface 7 is assigned the 
graph consisting of the trajectories n^^^ and n^^^ entering {x,i) and 
the trajectory n^f^ leaving this point (i.e., the trajectories n^J^ fiber 
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the surface F^'^). We also assume that the surface (stratified mani- 
fold) Fy = F(i) UF(2) UF(3) consists of points belonging to these graphs. 
Next, we assume that u{x,t) is a piecewise smooth vector field whose 
trajectories enter Fu. 
Definition 1.2. Let 

3 

u{x, t) = uo{x, + ^ H{t - 'i/ji)uii{x, t), 

i=l 

where -0 is the same function as before, and uo,uii G C(M"' x M^). The 
function p{x, t) determined by the relation 

3 

p{x, t) = R{x, i) + XI ^«(^)<^(* ~ 

i=l 

where R{x,t) G C^{W x R\) \ {[jvf^}, is called a generalized (5-shock 
wave type solution to (1.2) corresponding to the stratified manifold Fy 
if the integral identity 



1 / {RCt + (uR, VC) + aRQ dx dt 


3 



+E/ , ;,, c(^,t)dx=o (1.5) 



ei d_ 



«=i - - (i) 

holds for all test functions C(a;,i) G V{W x Rj^), C(^,0) = 0, = 

(^,V) + |V^,||. _ 

As above this relation implies the first equation from (1) outside Fy, 
equations of the type of the second equation in (1) on strata F(j-) and 
the Kirchhoff type relation on 7: 

(ei + 62)1- = 63!-. (1.6) 

Now we consider the case with codimF > 1. First, we note that the 
second integral in (1.2) can be written as 

We note that if the surface F is determined by the equation S{x, t) = 
rather than by the simpler equation {t = ip{x)} presented at the 
beginning of this section, then 

V _ vs _ s, 

|V5| \VS\~ |V5|2^^ 
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and, of course, the new vector field -^j^ = {Vn, V) + ^ remains tangent 
to r. 

Therefore, in this more general case, using this new vector Vn, we 
can again rewrite the integral identity from Definition 1.1 as 

/ / {RCt + {uR,VO+aRC)dxdt+ [ e({Vn,^) + -^)(:{x,t)dx = 0. 
Jo -/r" Jr ^ otJ 

(1-7) 

This form of integral identity can easily be generalized to the case 
in which F is a smooth surface in W"^^ of codimension > 1. 

In this case, instead oiVn, we can use a vector v that is transversal 
to r and such that the field (v, V) + ^ is tangent to T. We note that the 
vector V is uniquely determined by this condition, which can be treated 
as "the calculation of the velocity value on the discontinuity" from the 
viewpoint of [5] and [7]. 

Moreover, in this case, the expression for p does not contain the 
Heaviside function, and it is assumed that the trajectories of the field u 
are smooth, nonsingular outside F, and transversal to F at each point 
of F. In this case, the function p has the form 

p{x,t) = R{x,t) + e{x)5{V), 

where R G C^(M"+^\F), e G C^(F), and the function (5(F) is determined 

by 

(,5(F),C(x,t)) = ^Cc^, 

where oj is the Leray form on F. If F = {5*1 (x, t) = Ofl - • •n5fc(x, t) = 0}, 
k G [l,n], then oj is determined by the relation, see [6], p. 274, 

dtdxi dxn = dSi dSkCO. 

In this case, we assume that the functions are sufficiently smooth 
(for example, C^(M'^ x and their differentials on F are linearly 

independent. 

Moreover, we can assume that the inequality 

^_ ViS,,...,Sn) 

'D{t,Xi,. . . ,Xn-l) 

holds. This inequality is an analog of 5^ / at the beginning of this 
section and allows us to write u in the form 

OJ = J^^dxk dxn- 
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The integral identity, an analog of (1.7), has the form 
^ J^^ {RCt + {uR, VC) + aRC) dxdt + j^e(^{v,V) + ({x, t)oj = 0. 

Integrating the latter relation by parts, we obtain equations for deter- 
mining the functions e and R similarly to (1.4). 

Now we assume that the singular support of the velocity field is the 
stratified manifold UFj with smooth strata Tt of codimensions nj > 1. 

We also assume that the velocity field trajectories are transversal to 
U r and are entering trajectories. 

Then the general solution of Eq. (0.12) has the form 

p{x,t) = R{x,t) + Y.^i8{Ti), (1.8) 

where R{x,t) is a function smooth outside [jVi, are functions 
defined on the strata Fj, and the sum is taken over all strata Fj. 

The integral identities determining such a generalized solution have 
the form 

/ / (RCt + {uR, VC) + aRQ dxdt 
Jo Jw 

i -^^^ 

This implies that, outside IjFj, the function R satisfies the continuity 
equation 

Rt + {V,uR) + aR = 0, 

and, on the strata Tj for nj = 1, equations of the form (1.4) hold, 
which contain the values of R brought to F^ along the trajectories. For 
ni = n — k, k > 1, on the strata F;, we have the equations 

Q^eini + (V, vieini) = Fim, (1.10) 

where pi is the density of the measure with respect to the measure 
on F; which is left-invariant with respect to the field ^ + (f j , V) , and 
Fi is defined by the following construction. Denote a e-neighborhood of 
F; by F[ and denote its boundary by (?F|. Let us consider the integral 
appearing after integration by parts: 




(^^i,V) + ^ )C(x,t) 



(1.9) 
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where Uni is the normal component of velocity u on dFf , oof is the Leray 
measure on dVf, C, \s a. test function. Passing to the limit as e ^ we 
obtain 



It is well known that outside |J Fj the function R[x,t) can be calcu- 
lated using the famous Cauchy formula 



where po ia a constant along the trajectories of the field u outside IJ^i, 
\^§^\ is the jacobian of the mapping corresponding to the shift along 
the trajectories of u and the integral under exponent is calculating along 
the trajectories of the field u. 

This formula implies that the limit as e ^ of the above integral 
exists. 

We note that it follows from the above that the function R is deter- 
mined independently of the values of Vi on the strata under the condition 
that the field trajectories enter (jFj. 

In conclusion, we consider the case where the coefEcient a has a 
singular support on [jTi, i.e., 



where a = f{u) outside \J Fj. We note that such a choice of the definition 
of the term ap is not unique in this case. But, first, it is consistent with 
the common concept of measure solutions (see [3], [5]) and, second, it is 
of no importance for the construction of the solution outside |J Fj for 
the case in which the trajectories u enter IjFj. 
In this case, identity (1.9) takes the form 





(1.11) 



a = f{u). 



In this case, we set 





(1.12) 



and Eq. (1.7) can be rewritten in the form 
d 

-rrXeiiii) + {V ,vieini) + f{vi) = Fim. 



(1.13) 
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All the afore said gives the following statement. 



Theorem 1. Let that the following conditions be satisfied fort G [0,T], 



(1) IJ Fj is a stratifies manifold with smooth strata Fj/ 

(2) the trajectories of the field u are smooth outside [jFj, enter (jFj 
and do not intersect outside (JFj/ 

(3) equations (1-13) are solvable on the strata F,; 

(4) the Kirchhoff laws are satisfied on the intersections of strata Fj. 
Then there exist a general solution to the continuity equation (0.1) 

with a = f{u) in the sense of the integral identity (1.12). 

2 The Maslov tunnel asymptotics 

Wc recall that the asymptotic solutions of a general Cauchy problem for 
an equation with pure imaginary characteristics was first constructed by 
V. P. Maslov [8]. In the present paper, we consider only the following 
Cauchy problem 



-h^ + Prx,-h—]u = Q, u{x, t, h)\t=o = e-^°(-)/V°(a;), (2.1) 



where P{x,^^) is the (smooth) symbol of the Kolmogorov-Feller operator 
[9]) S'o > is a smooth function, ip^ G Cq°, h +0 is a small parameter 
characterizing the frequency and the amplitude of jumps of the Markov 
stochastic process having having transition probability given hy P{x,^). 
To be more precise, we can have in the mind the following form of 



where A{x) is positive definite smooth matrix and fj,{x, du) is a family 
of bounded measures smooth with respect to x. The symbol P{x, ^) can 
also depend on we will be more precise later on. 

Locally in t, an asymptotic solution of problem (2.1) can be con- 
structed according to the scheme of the WKB method, see [8]: the 
solution is constructed in the form 



T > 0: 



1 





oo 



U = 
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in the sense of asymptotic series. In this case, for the functions S{x,t) 
and ipo{x,t) we obtain the following problems: 

^ + =0, S{x,t)\t=o = So{x), (2.2) 



difo dS\ ^ \ d'^P d'^S , , 



dt 



ipo{x,t)\t=o = ^^{x), 



As is known, the solution of problem (2.2) is constructed using the 
solutions of the Hamiltonian system assumed to exist and to be smooth 

x = V^P{x,p), x|i=o = xo, (2.4) 

P = -'^xP{x,p), p|t=o = V5o(xo). 

This solution is smooth on the support of ipo{x,t) for all t such that 

the Jacobian Dx/Dxq ^ for xq G supp(^'^(x). We let g^^ denote the 

translation mapping along the trajectories of the Hamiltonian system 
(2.4). 

We recall that the plot 

K^ = {x = xq,p = VSq{xq)} 

is the initial Lagrangian manifold corresponding to Eq. (2.2), and AJ* = 

5(^Aq is the Lagrangian manifold corresponding to Eq. (2.2) at time t. 
Let TT : A" — )• be the projection of A" on M^, which is assumed to be 
proper. The point a G A" is said to be essential if 

S{a,t)= min S{^,t) 

/3e7r-i{a) 

and nonessential otherwise. Here S is the action on A" determined by 
the formula 

ft 



5(/3,t) = / pdx-Hdt, 
Jo 



where the integral is calculated along the trajectories of the system (2.4) 
the projection of whose origin is xq = (3. As is known 

S{x,t) = S{TT-^x,t) 

at regular points where the projection tt is bijective. 
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The global in time asymptotic solution of problem (2.1) is given by 
the Maslov tunnel canonical operator. 

To define this operator, following [8, 10] wc introduce the set of 
essential points Ij7i« <^ ^t- This set is closed because the projection vr 
is proper, i.e. that for all x the set of p such that {x,p) G A", 7r(a;,p) = x 
is finite. 

Suppose that the open domains Uj C A" form a locally finite covering 
of the set Ij7it- the set Uj consists of regular points, then we set 



(2.5) 



where 



Dxq 



Dx 



1/2 



where il^oj{x,t) being the solution of the equation 



^ + (P^(a;,V5,),VVo,; 



1 a^p 

- tr ^^^(x, VSj)'ipoj 



that exists and is smooth whenever 



dxd^ 



Dx/Dxq 



(2.6) 



7^ 0. The solution uj 



in the domain containing essential (nonregular) points (at which dir 
is degenerate) is given in the following way: the canonical change of 
variables is performed so that the nonregular points become regular, 
then we determine a fragment of the solution in new coordinates by 
formula (2.5) and return to the old variables, applying the "quantum" 
inverse canonical transformation to the solution obtained in the new 
coordinates. 

The Hamiltonian determining this canonical transformation has the 
form 

1 " 



where cJi, . . . C7„ = const > 0. 

The canonical transformation to the new variables is given by the 
translation by the time —1 along the trajectories of the Hamiltonian 
Ha- One can prove (see [8], [10]) that the family of sets a for which the 
change of variables takes a regular point into a nonregular is not empty. 

Next, the solution near the essential point is determined by the re- 
lation 

Ui = eTi^'"Uj, (2.7) 
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where Uj is given by formula (2.5) in the new variables and 




On the intersections of singular (containing singular points) and non- 
singular charts (without singular points), we must match Sj and ^qj- 
This can be done by applying the Laplace method to the integral whose 
kernel is a fundamental solution for the operator —h^ + Her. This in- 
tegral appears if we write the right-hand side of (2.7) in detail. In this 
case, since the solution is real, the Maslov index which is well-known 
[8] to appear in hyperbolic problems does not appear. The complete 
representation of the solution of problem (2.1) is obtained by summing 
functions of the type (2.5) and (2.7) over all the domains [/j, for more 
detail, see [8], [10]. 

The asymptotics thus constructed is justified, i.e., the proximity be- 
tween the exact and asymptotic solutions of the Cauchy problem (2.1) 
is proved [8, 9]. More precisely it is proved that at the points of the set 
7r([J 7it where the projection tt is bijective the following estimate holds: 

u{x, t, h) — Uj = 0{h) 

In the preceding case we noted that values of the solution of the 
continuity equation at nonregular points are independent of the values 

of the solution on the singularity support (of course, the inverse influence 
takes place) by the condition that the velocity field trajectories enters 
the singular support. 

In the case of the canonical operator construction briefly described 
above, the relation between the solutions at essential and nonessential 
point is also unilateral, namely, the essential points arc "bypassed" using 
(2.7), but the values of the functions V'oi contained in Uj on the singu- 
larity support do not determine the values at the regular points (but the 
converse is not true). 

Now we note that the function S{x, t) such that 

S{x,t)\u,=Sj{'K~\a),t) 

is globally determined and continuous at points of the domain 7r(|J 7^4) C 
M^. We denote this set by (jFj and assume that this is a stratifled 
manifold with smooth strata Ta of different codimensions. We note that, 
for example, if the inequality V{Si{x,t) — Sj{x,t)) 7^ holds while we 
pass from one branch AJ* n |J to another, then the set Tr{{Si — Sj) = 0} 
generates a smooth stratum of codimension 1. In the one-dimensional 
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case, all strata are points or curves on the (x, t)-plane (under the above 
assumptions about the singularities being discrete). 

Now we consider the equation for ipQj. We denote this function by p 
and then obtain 

^ + {V,up) + ap = 0, (2.8) 

where u{x, t) = VgP(x, VS) and a = — tr §^{x, VS). 
If the condition 

Hess^ P{x, > 

is satisfied, then it follows from the implicit function theorem that 
VS{x,t) = F{x,u{x,t)), where F{x,u) is a smooth function and 

a = f{x,u), 

where f{x, z) is again a smooth function. 

Let us return to the formula (1.8) and denote the regular (in the 
sense of distributions) part of p by preg- 

Thus, we can formulate the following theorem. 

Theorem 2. Suppose that the following conditions are satisfied for t G 

[o,r], T > 0.- 

(1) There exists a smooth solution of the Hamiltonian system (2.4). 

(2) The singularities of the velocity field u = V^P{x, VS) form a 
stratified manifold with smooth strata and Hess^ P{x,^) > 0. 

(3) There exists a generalized solution p of the Cauchy problem for 
Eq. (2.8) in the sense of the integral identity (1.10). 

Then at the points o/7r(|j7it) where the projection tt is bijective, the 
asymptotic solution of the Cauchy problem (2.1) has the form 

u = eyiY>{-S{x,t)/h){y/p^g + 0{h)). 

This theorem is a global in time analog of the corresponding Made- 
lung observation about local solutions of Schroedinger type equations. 

Now we demonstrate the connection between solutions of continu- 
ity and transport equations. It is easy to see that by construction a 
transport equation solution is equal to p^^g out of the points of singular 
support. More that one can define yfp globally in distributional sense 
and thus defined square root is equal to preg- To prove this statement 
it is sufficient to note that (5-shock type solution to continuity equation 
in the form (1.8) can be obtained [14, 15]) as a weak limit of a weak 
asymptotic solution of the form 

p{x, t) = p,eg + 5Z ('M^i)^ (2-9) 
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where x G i?", e is an axillary small parameter and is a regularization 
of 5-function of the form 

6, = £-("-"*)u;(5i/£). 

Here Si are smooth functions such that VS'j|5.=o / and the equation 
Si = defines the strata Fj, nj = dimFj and oj = uj{ri) belongs to the 
Schwartz space of test functions. Without loss of generality one can 
consider the case = in (x, t) half-plain for each fixed t. In this case 
the function Si can be chosen in the form S = x — (p{t) and our statement 
becomes 

l^^ojj^i ^'^'^"^ ei£-ia;((x - 0)/£)) - y/p^\ (p{x)dx = 0, (2.10) 



for arbitrary test function (p{x) from Cq 

After the change of variables rj = {x — (p)/e the (2.10) gets the equiv- 
alent form 

lim Vs j^^ (^^£preg + eiu;{ri) - ^/ep^ d-q = 0, (2.11) 

which is obviously true under the trivial inequality \/l + z < 1 + z/2. 
Indeed consider the domain in i?^ where < pe for some a, < a < 1. 
Then the mentioned inequality gives the estimate 0(e^/^) for the inte- 
gral under limit in (2.11). In the case uj°' > pe, because the function u is 
decreasing (belongs to the Schwartz space), the measure of correspond- 
ing domain can be estimated as {pe)~^/^ where N is arbitrary positive 
number. This completes the proof. 



3 Particular cases 

The theorem stated in the previous section requires that some assump- 
tions are satisfied. The most restrictive is the item 3 in the theorem 
above-that is the existence of the global generalized solution to conti- 
nuity equation. Under the above made assumptions it is possible to 
construct this solution using characteristics, but only in the case where 
the structure of singular support of u is not changing in time-all sections 
of the stratified manifold introduced above by planes t = const are dif- 
feomorphic. A more complicate situation arises, when the singularities 
of the velocity field change their structure. In this case the problem 
of the construction of a global in time generalized solution to the con- 
tinuity equation has not been solved yet. The obstacle is that in this 
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case usually one has no global in time expression for the velocity field u. 
In turn this does not allow to apply formula (1.11) to construct global 
solution to the continuity equation. In multy-dimensional case as far as 
we know there is only one result concerning to shock wave generation 
[14] which allows to construct global in time approximation of the shock 
wave formation process. But this is slight different from the construction 
that we needs here. In the one dimensional case the situation is better 
and we have all needed formulas. 

We begin with the spatially homogeneous case. Here the problem is 
equivalent to the one of constructing a formula for a global solution to 
conservation law equation 

Here P{-h-^,t) is the same operator as in (2.1) but assumed to 
be independent of x with the symbol P{S,,t) and v = dS/dx. The 
velocity field u in this case is P^{v, t). In [11] a construction of the global 
solution to the continuity equation where the velocity field is given by 
the solution of the equation (3.1) was given. Because the set of singular 
points is discrete by our assumptions, without loss of generality one can 
consider the case were only one point of singularity appears. Denote the 
corresponding (smooth) initial condition by uq, the instant where the 
singularity appears by t* and the point of singularity by x*. 

The first step of construction suggested in [11], [12] (see also [18]) 
is that we change uq in a small neighborhood of of origin Xq of the 
trajectory coming to x* when t = t*. We denote this new part of initial 
data as ui{xo) for xq G {xq — /3, Xq + /3), /3 — >■ and assume 

^> 0, e ^ 0. (3.2) 

We define the function ui = ui{xo,t) as a solution of implicit equa- 
tion 

Pl{ui,t) = -K{t)xo + b{t), (3.3) 

The latter equation is solvable under the condition Hess^ P(x, ^) > 0, 
formulated above. 

The functions K(t) and b{t) are defined from the condition of conti- 
nuity of the characteristics flow, i.e 

ui{xo - p,t) = Uo{Xq - P, t), Ui{Xq + p,t) = Uo{Xq + (3, t) 

It is easy to check that this choice of ui provides that the Jacobian 
|£)a;/I?a;o| is identically equal to for i = t* and xq G {xq — (3,Xq + 
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j3). Here we remove from usual topological concept of general position 
considering the situation of identical equality that can be destroyed by 
small perturbation. But this construction follows from the algebraic 
concept and allows to present the solution of (3.1) in the form of linear 
combination of Heaviside functions (see [11]). 

The second step of our construction of an approximation is a mod- 
ification of the definition of characteristics. We set 

x = {l-B)Pl{ui{xo,t),t) + Bc, xoe{x*o- P,x*Q + P) (3.4) 

and 



X = P^{uo,t), 
where xq does not belong to {xq — P,Xq + P), 

P{v{x{xl + /3, t),t)) - P{v{x{xl - /3, t),t)) 
v{xix*o + P,t),t)-vix{x*o-^,t),t) 

Initial data for (3.4) are the following: 

= xo + Ae, £ > 0. 

=0 

The function B in (3.4) has the form B = B{{t - t*)/e) and B(z) is 
smooth, monotone and increasing from to 1 for z G (—00,00). Simi- 
larly to [11], [12] one can prove that there exist an ^ = const such that 
the Jacobian |£)a;/-Dxo| calculated using the above introduced charac- 
teristics is not equal to zero, but it is of order 0{e) when t > t* + O{0) 
when Xq G {xq — P,Xq + /3). Using the velocity field generated by x we 
can construct global in time (smooth) solution of the continuity equa- 
tion in the form (1.11). After that, passing to the limit as e — > we will 
obtain the generalized solution of the continuity equation in the sense 
of definition from Sec.l just like it was done in [12]. 

Spatially inhomogeneous one dimensional case. 

We will follow the scheme introduced above. The case under con- 
sideration can be treated in the same way as the previous one with 
modifications. Firstly, we will assume that the symbol P = P(x,^) does 
not depend on t. In this case this assumption (which means that the 
mapping Qp is invertible)will be used to construct the insertion to initial 
data. In the previous case we did it using the implicit function theorem, 
see (3.3). 

Let Aq be a smooth nonsingular (w.r.t the projection tt) curve in 
the {x,p) space, which is a Lagrangian manifold corresponding to initial 
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data for our problem. We consider Lagrangian manifold A]* = gp Aq and 
assume that there is only one point singular with respect to projection 
on X-axis and its projection is x* . Let f3 be the same as above. Let 
us set tl = t* + (3. Because of the assumption that P^'^ is positive, we 
have that for t = the Lagrangian manifold Aj, has two parts which 
contain essential points and these parts form a shock wave type curve 
with the jump at the point xl where Si(,ft{xi,tl) = Srighti^i:ti)- Wc 
connect these parts by a vertical line and thus obtain a new Lagrangian 
manifold, which is piecewise smooth continuous curve with two angle 
points (ends of the vertical part, the distance between them of order 
/3). We denote this manifold by Aj. and apply the mapping (/p*^ for 
sufficiently small ti to this manifold. This mapping obviously exists and 
is a diffeomorpfism because our Hamiltonian P does not depend on t. We 
consider the obtained manifold gp^^ A|, as the new Lagrangian manifold 
corresponding to our problem for t = — ^ichanging the manifold A^*_^^ 

by 5p*^ A^* . As it was said above the latter manifold is piecewise smooth 
curve witii two angle points and all points of the curve outside of the 
part between these angle points are regular. Moreover there exist a 
sufficiently small ti such that the part of the curve between these angle 
points contains only regular points-these statements are the consequence 
of the positivity of P" , its stationarity and the possibility to choose ti 
small enough (and independent on s). 

Denote the projections of the mentioned above angle points on the 
manifold ^p*^ Aj* to the x-axis by ai < 02 and note that | ai — 02 1 is of 
order /3. 

Like in the previous example we introduce the new characteristics 
system 

X = (1 - B)P^{x{xo,t),p{xo,t)) + Be, (3.5) 
p = -{l- B)P^{x{xo,t),p{xo,t), xo e (01,02), 

and 



x = P^ix,p), (3.6) 

when Xo does not belong to (01,02). We have set 

P{v{x{a2,t),p{a2,t))) - P{x{ai,t),p{ai,t)) 



p{x{a2, t), t) - p{x{ai,t), t) 
Initial data for (3.5), (3.6) are the following: 



(3.7) 
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x\^^Q = xo + Ae, 
p\t^Q =Po{xo), 

where (xo,Po(a^o))=5p*^^j*- The expression in the right hand side of 
(3.7) is the direct analog of the well known Rankine-Hugoniot expression 
for the velocity of the shock propagation. In the case under consideration 
it is the velocity of the point x on x-axis, where Sieftix, t) = Srighti^, t). 

By the assumption we have only one singular point if we are consider- 
ing the family of manifolds Aj, t e [0,t*]. We also have by construction 
that the jacobian J = Dx/Dxq calculated using the solutions of the 
system (3.5) is not equal to zero. More precisely we have 

lim J = H{t* -t)Jo, 

where Jq is the Jacobian calculated using the solutions of (3.5) for i? = 
(Jo = when t = t* by construction) and 

J > H{t* - t)Jo + Ce, 

where C = const > 0. This statement directly follows from (3.5) if 
we take the properties of the function B into account. It means that 
the velocity field, generated by projections of the solution of the system 
(3.5), (3.6) on the x-ax;is has nonintersecting trajectories for e > 0. 
Thus we can use it to construct solutions of the continuity equation. It 
remains to note that just like in [12] it is easy co check the the limits of 
these solutions will satisfy to the integral identities introduced in Sec. 1 
as the definition of generalized solutions to continuity equation. 
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